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Abstract Normal copula with a correlation coefficient between —1 and 1 is tail independent 
and so it severely underestimates extreme probabilities. By letting the correlation coefficient in 
a normal copula depend on the sample size, Hiisler and Reiss (1989) showed that the tail can 
become asymptotically dependent. In this paper, we extend this result by deriving the limit of 
the normalized maximum of n independent observations, where the z-th observation follows from 
a normal copula with its correlation coefficient being either a parametric or a nonparametric 
function of i/n. Furthermore, both parametric and nonparametric inference for this unknown 
function are studied, which can be employed to test the condition in Hiisler and Reiss (1989). A 
simulation study and real data analysis are presented too. 
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1 Introduction 


Let {(X);, V)} 1 be independent and identically distributed random vectors with distribution function 
F(x, y), continuous marginals F\ and F 2 . The copula of F is defined as F(Ff (x), Fiff ( y )), where F~ denotes 
the inverse function of F). Assume the copula of F follows from a normal copula C(x, y ; p), where p £ [—1,1] 
is unknown. Hence the density of C(x,y;p) is 


c(x,y;p) 


1 (2p*-(x)*-(y) - p 2 (<S>~(x)) 2 - P 2 ($-(y)) 2 

yr^ 6XP V 2(1 -?) 


( 1 . 1 ) 


for p £ (—1,1), where $(x) is the standard normal distribution function. 


Normal copulas are one of most commonly used elliptical copulas, and elliptical copulas are popular 
in risk management due to their ease of simulation (see McNeil, Frey and Embrechts (2005)). Recently 
Cliannouf and L’Ecuyer (2012) used normal copulas to model arrive processes in a call center, Fung and 
Seneta (2011) showed that a bivariate normal copula is regularly varying, Meyer (2013) studied the properties 
of a bivariate normal copula, efficient estimation for bivariate normal copula models was studied by Klaassen 
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and Wellner (1997). Although normal copulas are easy to use and have some attractive properties, a serious 
drawback of using a normal copula is the so-called tail asymptotic independence (see Sibuya (I960)), which 
under-estimates extreme probabilities in risk management. 

To overcome the shortcoming of the tail asymptotic independence of a normal copula, Frick and Reiss 
(2013) assumed that p = p(n) satisfies the so-called Hiisler-Reiss condition 


(1 - p) log n 

(cf. Hiisler and Reiss (1989)) and proved that 


A € [0, oo] as 


^n(maxi<j<„ Fi(Xi) - 1) < x, n(maxi<j<„ F 2 (Yi) - 1) < yj 

($(Va + ^0M)x + $(Va- 


exp 


iog(y/a:) 

2%/X 


( 1 . 2 ) 


(1.3) 


for x < 0 and y < 0 as n —> oo. This is the copula version of the limit in Hiisler and Reiss (1989) for the 
normalized maxima of n independent random vectors with a bivariate normal distribution and its correlation 


coefficient satisfying (1.2). Obviously, a bivariate random vector with the above limiting distribution is 
dependent when A < oo. Extending the results in Hiisler and Reiss (1989) to elliptical triangular arrays is 
given in Hashorva (2005, 2006). 

Since the above p depends on the sample size n, one may call it dynamic normal copula. Recently dynamic 
copulas are receiving some attention in modeling financial time series; see Benth and Kettler (2011), Mendes 
and de Melo (2010), Guegan and Zhang (2010), and Van den Goorbergh, Genest and Werker (2005). 


In this paper, we further study the convergence in (1.3) by allowing p to depend on both i and n. That 


is, we do not assume that (V, : ,Y))'.s are identically distributed. Motivated by (1.2), an obvious extension is 
to assume that (1 — p) logn is a function of i and n. As in nonparametric regression models, we assume that 
(1 — p) log?r is a smoothing nonparametric or parametric function of i/n so that we can employ well-developed 
local polynomial techniques to estimate this function and to test whether this function is a constant, which 
gives a way to verify the condition imposed by Hiisler and Reiss (1989) and Frick and Reiss (2013), and 
indicates the observations have the same distribution. More specifically we assume that {(Xj,l))}" =1 is 
a sequence of independent random vectors and the copula of (X, ,l)) is a normal copula with correlation 
coefficient p = 1 — m(i/n)/ log n for an unknown smooth function m( x). After deriving the convergence 
for the normalized maxima of the copulas of (Xj,Y)) / s, we propose both parametric and nonparametric 
estimation for m(x), which are based on either Kendall’s tau or correlation coefficient. We also derive the 
asymptotic limits of the proposed estimators, which turn out to be quite nonstandard with an unusual rate 
of convergence. The proposed estimators can be used to determine tail dependence, which is of importance 
in predicting co-movement in financial markets; see McNeil, Frey and Embrechts (2005). 

We organize this paper as follows. Section 2 presents the main results and statistical inference procedures. 
A simulation study is given in Section 3. Section 4 reports some empirical data analyses. All the proofs are 
given in Section 5. 
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2 Methodology 


Throughout, suppose {(Xj,y )}" =1 are independent random vectors, X's have the same continuous dis¬ 
tribution function F± and Y- s have the same continuous distribution function F 2 . Assume the copula of 
(Xi , Yi) is the normal copula C(x, y; pi) with density given by & 


2.1 Convergence of maxima and tail coefficient 


As motivated in the introduction, we extend the result (1.3) by assuming 


Pi = 1 — m(i/n)/\ogn for some nonnegative function m(s), 


which includes condition ( 1 . 2 ) as a special case. 


Theorem 1. Under condition (2.1), 


i) */maxi<j<„ m{i/n) —> 0 , then for any x < 0 and y < 0 

lim P ( n( max FAXA — 1 ) < x, n( max F 2 (Yi) — 1 ) < y ) = exp ( min(x, y) ); 

n—foo \ l<i<n l<i<n / \ / 

ii) if m(i/n ) — > oo, then for any x < 0 and y < 0 

lim P (n( max Fi(Xj) — 1) < x,n( max F 2 (Yi ) — 1) < y) = exp(x + y); 

n—>oo \ l<i<n l<i<n / 

in) if m(s) is a continuous positive function on [ 0 ,1], then for any x < 0 and y < 0 

lim„^ oc P ^n(maxi<j<„ Fi(Xj) - 1) < x, n(maxi<j<„ F 2 (Yi) - 1) < yj 

= exp(x/ 0 1 $(^(i) + g«) ds + yf^^V^+T^) ds^ 

=: G(x,y). 

Furthermore the tail dependence function l(x,y) = lim^o i _1 {l — G(tx,ty)} equals 

-x ( $ ds -V I $ (VMs) + ds 

Jo \ 2 ym(s) j Jo \ 2 yro(s )) 


( 2 . 1 ) 


for x < 0 and y < 0, and the tail coefficient is A = /(—1, —1) = 2 f* $(i/m(s)) ds. 


2.2 Parametric inference 

Here we consider statistical inference for fitting a parametric form to the unknown function m(s). First, 
we consider the family m(s) = a + (3s 1 , where a > 0, (3 ^ 0 ,7 > 0. Note that when /3 = 0, 7 can not be 
identified. Also when 7 = 0, a and (3 can not be distinguished. 

It follows from Theorem 5.36 of McNeil, Frey and Embrechts (2005) that 

E [sgn({Ui - Ui)(Vi - V*))) = ^ arcsin(pj), (2-2) 
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where (Ui, V.) is an independent copy of ( Ui , Vi), and 


1 


1 


1 


E ( (Ui - ^)(Vi - 7) ) = A- arcsin(^). 


Also we have 


(2.3) 

(2.4) 


E($-(F 1 (X i ))$-(F 2 (y i ))) = Pi- 
Therefore, one can employ the standard least squares estimate based on one of the above equations. 

Since ( Ui,Vi)'s are not identically distributed, we do not have an independent copy of (t/j,V)), which 
prevents us from using (2.2). Hence we propose to use either (2.3) or (2.41 to construct the least squares 
estimator, which results in 

(a 


n / i | 

>£,7) = arg min ^ (ApQ) - -)(F 2 (Yi) - -) - 
(a,^, 7 ) ^ \ Z Z 


1 . 1 - (a + /3(z/n) 7 )/logn 

— arcsmf-) 

2 tt v 2 ' 


or 


n , 

(d*,/3*,r) = arg min V $~(F 1 (X i ))$-(F 2 (y i )) - 1 + 

(a,/3, 7 ) ^ \ 


/3(i/n)' 


logn 


where F^x) = ^ H X i < x ) and F 2 (y) = ^ =1 /(y 4 < y). Alternatively we define (a,j3,j) to 

be the solution to the following score equations 


Zni(a»/3,7) 
ln2(oi, (3, 7) 
ln3{OL,P,i) 


= E"=i ((A(A,) - \){F 2 {Yi) - |) - £ arcsin(f)) = 0, 

= Er=i ((A(AO - i)(F 2 (y,) - i) - A arcsin( f ))(i)7 = 0 , 


2tt ' 
1 


(2.5) 


= E"=1 (fiW - §)( W) - |) - Aarcsin(f) (^)'nog(i) = 0 


and (d*,/3*,7*) to be the solution to the following score equations 


Ci(<T £, 7 ) 


= E?=i ($-(A(Ai))$-(F 2 (li)) -1 + ° +w/n) A - 


= E"=i f$-(A(A i )) 4 >-(y 2 (y l )) -1 + 
= E?=i (^-(F^Xi^-iF^Yi)) -1 + 


log n 

x -\-(3 (i / n ) 1 
log n 

o:+/3(z/n)' 7 
log n 


= o, 

(A 7 = o, 

(irio g (i) = o. 


( 2 . 6 ) 


Note that we skip the term of A- arcsin(p,/2) in (2.5), which goes to a constant uniformly in i since pt —► 1 
uniformly in i. 


The following theorems give the asymptotic normality of the proposed estimators. 

Theorem 2. Suppose \2.1\) holds with m(s) = a + /3s 7 for some a > 0, /3 7 ^ 0 ,7 > 0. Then we have 


N( 0,E) 


/ V" 

0 

° 1 


^ a — a\ 

(log n) 3 / 4 


0 

y/n 
log n 

0 

A 


V 0 

0 

Y™ 

log n/ 


VT-7 / 


and 


where 


( yjr> 


V log n 


(■a — a), 


logn 


logn 


n 0-P), - 7 )) A /v(o, a- 1 e 0 (A~ 1 ) 


— 1\T 


(2.7) 


( 2 . 8 ) 



An 

0 

0 ^ 


(0 

0 

0 ^ 

S = 

0 

022 

023 

j E 0 = 

0 

022 

023 


1° 

023 

0"33 / 


1° 

023 

0"33 / 


4 

























an — v2 
cr 22 = 


y/o. + /3s 7 ds'j ^ fg(u - |) 2 </>(<1> (u))duj, 


&23 = 


180(1+27) 180(l+7) 2 ’ 

./j.. . +A 1 

180(l+27) 2 


0"33 — 


90(l+2 7 ) 3 180(l+7) 4 ’ 




1 

1+7 


1 _ §_ 

1+7 (1+7) 

1 0 

l+2 7 (l+2 7 ) 2 

2/3 


180(l+7) 3 

\ 


1 1 2/3 

\ (1+7) 2 (1+27) 2 (1+27) 3 / 


A = # 

67T 


/ 


1 

1+7 


1 

1+7 

1 

I+27 
1 


\ 


(i+W 

0 

(1+27) 2 

2/3 

(I+27) 3 / 


V (1+7) 2 (1+27)2 

Theorem 3. Suppose \2.1\) holds with m(s) = a + /3s 7 for some a > 0, /3 7 / 0, 7 > 0. TTien we have 


0 

y/n 
log n 

0 


y/n 
log n , 


\ 

a* — a\ 

A* 


/ 

+* - 7 ) 


and 


where 


( V™ 


\logn 


(a* - a), 


log n 


n - 1 ) 


logn 


7V(0, S*) 


N(0, (A*) _ 1 Sq(A * 2 ) _1 ) 


(2.9) 


( 2 . 10 ) 



/(T ii 

0 

0 N 


(0 

0 

0 > 

E* = 

0 

_* 

°22 

_* 

°23 

y* _ 

5 ^0 — 

0 

_* 

°22 

_* 

°23 


1° 

_* 

°23 

°33/ 


1° 

_* 

°23 

a 33/ 


'11 


2 33 


= 4(a + T ^) 2 + 2/3 2 ( 

4 2 


I+27 


r), 


( 1+7 ) 2 + 
2 


2 22 


1+27 (1 +7) 2 ’ 


+ 


(l+2 7 ) 3 (1+7) 4 ’ 23 (1+27) 2 (1+7) 3 ’ 

A* = 27 / 371 A and A* = 2 7/371 A, where A and A are given in Theorem^ 

Remark 1. Since 0-22 < °22/(127i 2 ) and <733 < er3 3 /(127i 2 ) ; /3 and 7 /iare a smaller asymptotic variance 
than (3* and 7 *, respectively, while the comparison for the asymptotic variances of a and a* is unclear since 
both — a) and — a) converge in distribution to zero. On the other hand, if one is interested 

in estimating A(a, f3,'y) T , then the estimator for the first element based on (d*,/ 3 *, 7 *) T has a faster rate of 
convergence than the corresponding estimator based on (d,/3, 7 ) T , but the estimators for the second and third 
elements based on (a*,/ 3 *, 7 *) T /iare a larger asymptotic variance than those based on (a,/3, 7 ). In spite 
of these theoretical comparisons, the simulation study below does prefer the estimation procedure based on 
equation (2.3) when the mean squared error is concerned. For testing (a,/3, 7 + = (ao, A), 7 o) T > one should 
employ the well-known Hotelling T 2 test statistic based on either (2.1) or (2.9) because the limit in both (2.8) 


and (2.10) is degenerate. 


Another interesting parametric form for m(s) is polynomial. Here we consider m(s) = a + /3s. In this 
case, when (0 = 0, m(s) becomes constant, which means that the observations (Xi,Yi), - ■ ■ ,(X n ,Y n ) are 
independent and identically distributed random vectors. 


Theorem 4. Suppose j[2.1\ ) holds with m(s) = a + (Is for some a > 0, (3 £ R. Then we have 

/3 


V™ , fi (3^ 

(logn ) 3 / 4 a 2 “ ^ 


n .a 

2 n logn 2 + 3 
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where E = 12n 2 (aij) is a symmetric matrix with 


(Tn = 


2y/2 

W 


((a + /3) 3/2 -a 3/2 )^ (u- *) 2 <K$ (w))d«. ^22 = £12 = 0. 


Theorem 5. Suppose \2.1\) holds with m(s) = a + f3s for some a > 0,/3 € R. Then we have 


rt~* , P E Vn A* , P* 

Vn{a + ~ a - I tt 

2 2 log n 2 3 


where E* = (a* ) is a symmetric matrix with 


IB 2 

ah = 4 a 2 + 4af3 + —— 
o 


a (3 


N( 0, E*), 


d 2*2 =g, *12 = 0 . 


Remark 2. When (2.1) holds with m(s) = a, we can show the rate of convergence for a* is faster than the 


rate of convergence for a. That is, the estimator based on (2-4) is preferred to that based on (2.3). However, 


the simulation study below prefers the estimation procedure based on equation (2.3) when the mean squared 
error is used as a criterion. 


2.3 Nonparametric inference 


First we use (2.3) to estimate the smooth function Q(s) = ^ arcsin( 

Especially we consider the local linear estimator Q(s) defined as 

{Q{s),b) = argmin^ f(ApQ) - \){F 2 {Yi) - - a - b{s - i/n) \ k( S ^ U ), 

a ’ i—1 ' ' 

where k is a kernel function and h = h(n) —> 0 is a bandwidth. That is, 

e;=i ^(A(^)-i)(F 2 «)-i) 


nonparametrically. 


Q(s) = 


E n 

i=i ’ 


where Wj = k( s J j n )[s na - (s - j/n)s n> i], s n< i = E"=i H * ~ j A) 1 ■ We refer to Fan and Gijbels [3 

for details. Therefore we can estimate m(s) non parametrically by 

m(s) = ^1 — 2sin(27rQ(s))^ logn. 

Theorem 6. Assume k(s) is symmetric with support [—1,1]. For a given s € (0,1), assume m"(t) is 
continuous at s, h = h(n) —► 0 and \((^ -> A as n ^ oo. Then as n —> oo we have 

Vnh , 


log?r 


(fh(s) — m(s)^j A N J t 2 k(t)dt , E J k 2 (t)dt^J . 


Second we use (2.4) to estimate the smooth function m{s) nonparametrically by considering the local 
linear estimator 


i—l 


(m*(s), 6) = argmin^ (ApG))$ {F 2 {Yi)) - 1 + ^—+ j—(s - i/n)) k( S 


i.e., 


m*(s) = —- 


log n log n 


E;=i^(^-(A(^))$-(F 2 (E)) - 1) logn 


, s — i/n , 
~h 


E TL 

j— i W 1 
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Theorem 7. Assume k(s) is symmetric with support [—1,1]. For a given s £ (0,1), assume m"(t) is 
continuous at s, h = h(n) —>■ 0 and h l (f^‘ —> A asn-> oo. Then as n oo we have 



logn 


(m*(s) 


m(s)) -4- N ^Am"(s) J t 2 k(t)dt, 



Remark 3. It follows from Theorems [d| and[?| that both rh*(s) and m(s) have the same asymptotic bias, but 
m(s) has a smaller asymptotic variance than m*(s). Hence, unlike parametric estimation, nonparametric 
estimation based on (2.3) is always preferred. 


Remark 4. By minimizing the asymptotic mean squared error, the optimal choices of h for rh(s) and rh*(s) 
are 

7r 2 k 2 (t) dt 


and 


hn 


: _ 
"0 — 


, 2x1/5 

log n 


1/5 


15(m"(s) f^tfkft) dt) 2 / 


log 2 n\ ( 2 f_ 1 k 2 (t) dt 


1/5 


_ (m"(s) f_ 1 t 2 k(t) dt) 2 / 


respectively, which are different from the standard optimal order 5 in the bandwidth choice of nonpara¬ 
metric regression estimation and nonparametric density estimation. Data driven method for choosing the 
above ho and can be obtained via estimating m"(s). A future research is to investigate the possibility of 
using cross-validation method to choose the optimal bandwidth. 


Remark 5. It is straightforward to construct both parametric and nonparametric estimation for the tail 
dependence function and the tail coefficient given in Theorem [7] and to derive the corresponding asymptotic 
results by using Theorems [1][7[ 


3 Simulation 

In this section we examine the finite sample performance of the proposed estimators by drawing inde¬ 
pendent (Xi,Y\), ■ ■ ■ , (X n ,Y n ) with ( Xi,Yi) following the normal copula with correlation coefficient p = 
1 — m(i/n)/logn. We consider n = 300 or 1000 or 3000, and repeat 1000 times. 

First we consider m(s) = a with a = 1 or 10, and calculate the average, sample variance and mean 
squared error for both a and a*. Table 1 below shows that a* has a smaller variance than d, which confirms 
the argument mentioned in Remark 2 that estimator a* has a faster rate of convergence than a. We also 
observe from Table 1 that i) a* has a larger bias and a larger mean squared error than a except the case 
of a = 10 and n = 3000; ii) the variance and mean squared error of both a and a* become larger when a 
increases; iii) the accuracy for both estimators improves as n becomes larger. In conclusion, a has an overall 
better finite sample behavior in terms of mean squared error than a* although its asymptotic variance is 
larger theoretically and empirically. 

Next we consider the case of m(s) = a + (3s. In Table 2 we report the average, sample variance and 
mean squared error for estimators (d,/3), (d*,/3*), (a + f + (|) and (a* + \ As we see, 

estimators (a, f3) have a smaller variance than (d*, / 3*), but d* + 4j- has a smaller variance than a+ which 
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Table 1: Estimators for the case of m(s) = a. 



a = 1 

a = 10 

a = 1 

a = 10 

a = 1 

a = 10 


n = 300 

n = 300 

n = 1000 

n = 1000 

n = 3000 

n = 3000 

E(d) 

1.0365 

9.9660 

1.0145 

9.9836 

1.0065 

9.9976 

V(a) 

0.0144 

0.0249 

0.0045 

0.0384 

0.0016 

0.0218 

MSE(d) 

0.0157 

0.0203 

0.0047 

0.0387 

0.0016 

0.0218 

E(d*) 

1.1690 

9.8440 

1.0788 

9.9476 

1.0368 

9.9914 

V(cU) 

0.0109 

0.0191 

0.0034 

0.0322 

0.0012 

0.0193 

MSE(d*) 

0.0395 

0.0434 

0.0096 

0.0349 

0.0026 

0.0194 


is supported by Theorems [ 4 ] and [5] that a* + has a faster rate of convergence than a + |. As n becomes 
larger, the accuracy of all estimators improves. Since a and /3 have a smaller mean squared error than a* 
and /3*, respectively, we prefer the estimation procedure based on equation (2.3) to that based on equation 


(2.4). 


Finally we consider the case of m(s ) = a + /3s 7 . Given results in Tables 1 and 2, we only consider the 


estimators derived from equation (2.3) with the large sample size n = 3000. Table 3 shows that all estimators 
have a rather large variance for 7 = 1 , and the variance of 7 is still quite big even when 7 = 0.5, which 
means estimating the shape parameter 7 is very challenging as usually. 


4 Data Analysis 

In this section we apply the proposed nonparametric estimators to two real data sets: Danish fire loss 
and log-returns of exchange rates; see Figure 1. 

This first data set is the nonzero losses to building and content in the Danish fire insurance claims, which 
comprises 2167 fire losses over the period 1980 to 1990. The second data set is the log-returns of the exchange 
rates between Euro and US dollar and those between British pound and US dollar from January 3, 2000 till 
December 19, 2007. 

We calculate both m(s) and m*(s) for s = 0.1,0.11,0.12, • • • ,0.9 by using Epanechnikov kernel k(x) = 
|(1 — x 2 )I(\x\ < 1) and the bandwidth h = d{log 2 (n)/n } 1 / 5 with d = 0.2,0.3, 0.4,0.5. From Figures 2 and 
3, we observe that m(s) and m*(s) have a quite similar pattern for the second data set, but seem having a 
different pattern for the first data set when a large bandwidth is employed. To further investigate this issue, 
we plot the difference of rh(s) — rh*(s) in Figure 4 for the above h with d = 0.2, 0.3, 0.4,0.5, which indeed 
shows the differences for d = 0.4 are quite similar to those for d = 0.5. Nevertheless, Remark 3 says that 
one should prefer m(s) to m*(s). The non-constant m(s) function indicates observations are not identically 
distributed. 









Table 2: Estimators for the case of m(s) = a + /3s with a = 1. 



/3 = 1 

/3 = 0 

0 = 1 

13 = 0 

0 = 1 

0 = 0 


n = 300 

n = 300 

n = 1000 

n = 1000 

n = 3000 

n = 3000 

E(d) 

1.0289 

1.0270 

1.0350 

1.0266 

1.0019 

1.0000 

V(d) 

0.2901 

0.3230 

0.1245 

0.1327 

0.0453 

0.0528 

MSE(d) 

0.2909 

0.3237 

0.1257 

0.1334 

0.0453 

0.0528 

E(/3) 

1.0345 

0.0486 

0.9612 

-0.0240 

1.0022 

0.0111 

V(/3) 

1.1597 

1.2678 

0.5080 

0.5095 

0.1833 

0.2082 

MSE(0) 

1.1609 

1.2702 

0.5095 

0.5101 

0.1833 

0.2083 

E(d+f) 

1.5461 

1.0513 

1.5030 

1.0146 

1.5030 

1.0055 

v(d +1)_ 

0.0270 

0.0150 

0.0097 

0.0047 

0.0032 

0.0015 

MSE(d+ |) 

0.0291 

0.0176 

0.0097 

0.0049 

0.0041 

0.0015 

E(f+ |) 

0.8593 

0.5297 

0.8379 

0.5053 

0.8350 

0.5037 

V(f + f)_ 

0.0170 

0.0131 

0.0070 

0.0047 

0.0024 

0.0019 

MSE(f+ |) 

0.0177 

0.0140 

0.0070 

0.0047 

0.0024 

0.0019 

E(d*) 

1.1654 

1.1557 

1.1155 

1.0880 

1.0349 

1.0292 

V(«*) 

0.4281 

0.4482 

0.2303 

0.2246 

0.0934 

0.1121 

MSE(d*) 

0.4555 

0.4724 

0.2436 

0.2323 

0.0946 

0.1130 

E(/3*) 

0.9802 

0.0338 

0.9188 

-0.0177 

0.9931 

0.0147 

v(4*) 

1.6477 

1.7563 

0.9138 

0.8838 

0.3686 

0.4504 

MSE(/3*) 

1.6481 

1.7575 

0.9204 

0.8841 

0.3686 

0.4506 

E(d* + £) 

1.6555 

1.1726 

1.5749 

1.0792 

1.5315 

1.0365 

V(d* + 

0.0221 

0.0112 

0.0076 

0.0037 

0.0025 

0.0012 

MSE(d* + Si) 

0.0463 

0.0410 

0.0132 

0.0100 

0.0035 

0.0025 

E(f + t) 

0.9094 

0.5891 

0.8640 

0.5381 

0.8485 

0.5195 

V(^ + ?) 

0.0174 

0.0152 

0.0087 

0.0071 

0.0033 

0.0036 

MSE(^ + 

0.0232 

0.0231 

0.0096 

0.0086 

0.0035 

0.0040 


9 



Table 3: Estimators for the case of m(s) = a + fis 1 with a = f3 = 1. 


E(d) 

V(d) 

MSE(d) 

E(/3) 

V(/3) 

MSE(/3) 

e(tO 

V(7) 

MSE(7) 

7 = 0.5 0.8631 

0.2050 

0.2237 

1.2268 

0.2840 

0.3354 

0.9787 

8.0693 

8.2985 

7 = 1 0.9964 

15.8532 

15.8532 

1.1412 

16.2379 

16.2578 

1.7859 

11.1177 

11.7353 


Danish fire 









s J 



Figure 1: Left panel: Danish fire loss with 2167 fire losses over the period 1980 to 1990. Right panel: 
log-returns of exchange rates between Euro and US dollar and those between British pound and US dollar 
from January 3, 2000 till December 19, 2007. 






Figure 2: Danish fire losses. Solid line and dotted line represent m(s) and m*(s), respectively. Bandwidth 
h = d{log 2 (n)/n} 1 / 5 with d = 0.2, 0.3, 0.4, 0.5 is employed in the upper left, upper right, lower left, lower- 
right panels, respectively. 

5 Proofs 

Proof of Theorem 1. We focus on the proof of case iii) since the other two cases can be verified easily. 
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Figure 3: Exchange rates. Solid line and dotted line represent m(s) and rh*(s), respectively. Bandwidth 
h = d{log 2 (n)/n} 1 / 5 with d = 0.2, 0.3, 0.4, 0.5 is employed in the upper left, upper right, lower left, lower 
right panels, respectively. 



0.2 0.4 0.6 O.S 


Figure 4: Danish fire losses. Differences of rh(s) — m*(s) are plotted for bandwidth h = d{log 2 (n)/?r} 1 / 5 
with d = 0.2, 0.3,0.4,0.5. 


For any e > 0 such that e < —x, write 


1 - P ^-(F^Xi)) < $-(1 + f ),$-(F 2 (y i )) < $-(1 + v n )) 

(FiiXi)) > $-(i + 2 ),$-(F 2 (F i )) > $-(i + 1)) 
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y 

n 

n 

X 

y 

n 

n 

X 

y 

n 
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/ ■oo 
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<£ - ( l -\- x / r . 




0-K 3a 7fir !£ )) 


$ (l+y/ra)-pjS N 

a/w! 

(l+i//ra) —(l+s/ra) ' 

\A-p? 

-1 f* (1+y/ri,)—(1+s/ri) \ , _1 f-e * / 4>“ (l+y/n)- Pi <S>~ (1+s/n) \ > 


d$(s) 
ds 


(5.1) 














For fixed x < 0 and y < 0, we have 


and 


K-n L / \ /hi- lo g (~y) loglogn + log(47r) log log n 

$ (l + 2/ /n) = y21^-- 7 _- 2v /2log^ + ° ( 7iog^ ) 

K-n , / >, ah - lo g(- s ) log log n + log(4?r) loglogn 

* (1 + s/n) = v ^"-7S5I5- 2V21ogn + ° ( vispr ) 


uniformly in s £ [x, — e], which implies that 


$ (1+y/n)-pj$ (1+s/ra) 

\A-pf 

_ V2 log riVl-pi _ log(-y) j _ Pi log(-s) _ 

Vl+Pi \/2 log ny/l—piy/l + Pi y/2 log Uyjl — pi y /1 +pi 2y/2\ogU </1 + pi ' ^ ^ \/log Ti 

\/ 2m(i/n )_ log(—p) _. (1 — m(i/n)/ log ra) log( — s) , /loglogn 


log log n+log(47r) <yi-pi . / log log n \ 
' \/Tnfrri. ' 


___(_ ^/- logiogn -v 

yj2—m(i/n)/ log n yj2m[i/n)2 — m(i/n) / log n 2m(i / n)yj2—m{i jn) / log n V^g n 

uniformly for s £ [x, —e], where x < 0 and y < 0 are fixed and e £ (0, —x) is any given constant. 
Since m(s ) is a continuous positive function, it follows from ( |5.2[ ) that 

„-i f x <&( — ~( 1 +v/ n )-Pi^ ~(i+a/n) \ j 

J - £ V v/wf / 

f x <&(\/m{i/n) - l °si~y) _|— i°g(-s) A (i -)_ o(l)) 

^J_ e ) 2v /m(i/n) ^ 2 y/m(i/n)) ) K K>> 


= (n 1 x4>f x /m(i/n) + lo ^ a / y ) ) + n 1 e$f yjm{i/n) + los / ) 

y V 2 y/m(i/n)J V 2y m(i/n) / 

—nr 1 — . 1 E ij>(Wm(i/n) - lo g(~n) _|— lQ g(~ s ) A n _|_ o(l)) 

2^/m(l/n) ^ e 2^/ m(i/n ) 2 y^m(i/n)J J 

= ( / n~ 1 x4> r m(i/n) + + n~ 1 e$(v / w(i/n) + 

y \ V 2i Jm(i/n) J \ V 2 ym(i/n) / 


+n 


2 y/m(i/n) ^ lo S 


log(-x) 


< b(y/m{i/n) 


log (~ 2 /) 

2^/ m(i/n 


2 / V/ i ) eSds ) (! + °( 1 )) 

2 y/m{%/n)J J 


= (n 1 x4>f Jm{i/n) + lo ^ 3 '/^ ) + n 1 e$f Wm(i/n) + lo ^ e,/y ^ ) 

\ V V W ' 2 y/m(i/n)J V v W ' 2 y/m(i/n) J 

- \/ rn ( i / n j) + ™ _1 y$(^== ~ VW*/™))^ (! +°( 1 )) 

= fn~ 1 x4>f -Jm{i/n) + + n~ 1 e4>f Jm{i/n) + 

y V 2-y m(i/n) J \ v 2 y/m(i/n)J 

-n-^y + rr^^^L + vWi/n)) “ V^AO)) (1 + «(!)). 

where 0(s) = ^'(s). Hence, it follows from (5.1) and (5.3) that 

lim^olmWoo £?=i (l-P($-(Fi(Xi)) < $~(1+ ^),$-(F 2 (^;)) < $"(1 + 


(5.2) 


(5.3) 


for any x < 0 and y < 0, which implies that 

lim^oo P ^(maxi<K n Fi(Xi) - 1) < x, n(maxi<j<„ F 2 (Fj) - 1) < 

= lim^oo n” =1 p (^"(^(Xi)) < $-(1 + z),*-{F 2 (Yi)) < $-(1 + f)) 

= exp flim^ Er=i logP ($-(Fi(^)) < Ml + f ),$-(F 2 (F0) < $"(1 + £))) 

= exp (- iim MOO E”=i (i - p (VCFiPQ) < *“(i + f), $-(F 2 (y,)) < $“(1 + »)))) 
= exp(x/>( v ^+^g)rf S + y/ 0 1 $( v 4^) + ^«)d S ) 
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for all x < 0 and y < 0. The rest for computing the tail dependence function and tail coefficient is 
straightforward. Q 


Proof of Theorem [§ Put Ui = F 1 (X i ), V z = F 2 (Yi ), U n (u) = ^ E"=i V u i < u )> V n (v) = ;ppr E"=i V v i < 
u), = (Ui - |) (Vi - 5 ) and = (u n (Ui) - 5 j (v n O'*) - 5 ) for * = 1 , • ■ • , n. Then 


- I ) ( UYf) - \ 


It is also known that 


sup 

0<w<l 


y/n{U n (u) - u} 


5 (1 — u) 6 


= O p ( 1 ) and sup 


0<u<l 


= {Zi}? = i- 


Vn{V n (v) - u} 


5 (1 — v) s 


— O p ( 1); 


see Inequality 1 in Page 134 of Shorack and Wellner [171(1986). 

Put 

1 n 

h = ~r Y,{u n (Ui) - UiXKiVi) - Vi), 

„ /on Z- J 


i= 1 


h = 


h = 


and 


p+tk/s Ei .1 E".! (mu, <u,)-u,)(v,-\) 

-fo So (HU, <u)~ u)(v - l)c(u,v;Pi) dudv), 

isdiTs Er.i E”„ (m < - v,m - 1) 

-fofim <v)-v)(u- \)c(u,V,p , 

Zi = (^T) E”=i fo JoWi <u)~ u)(v - \)c(u, v; Pj ) dudv 
+ (^T) E”=i fo 1 JoiW < v )~ v )( u - |) c ( u , V Pj) dudv 
+ (e ~ 5 “ arcsin(^)^ 

— Zf 1 Zi 2 “t - Z1 3. 


Therefore 


1 n 1 1 n 

arCsin (f )) =/l+/2+/3 + ^=5>. 

v i—1 v 2—1 


It follows from (5.5) that 


Ii — O p ( /—). 


Direct calculations show that E/f = O(f) and Eif = O(^), which imply that 


By (5.41, (5.6)—(5.8), we have 


Using 


h=O p (Z=) and I 3 = O p (-^=). 
\/Ti \/n 


1 1 ' ^ ~ _ 
—lni(a,f3,^) = —=y2Zi + O p (l/dJn). 

Vn Vnf-{ 


du 


d_ 

dv 


C(u,v;pi) = ~Ci(u,v,pi) 

C(u,v,pi) = <&( $ ft) ) := C 2 (u,v;pi), 


(5.4) 


(5.5) 


(5.6) 

(5.7) 

(5.8) 

(5.9) 
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we have 


and 


f 0 ( v ~ l)c(u,v,pi)dv 
= fn( v ~ \)Ci{u,dv\pi) 


= \Ci(u , 1; pi) + |Ci(it, 0;pi) - f 0 C 1 (u,v;p l ) dv 


f U (f){ £ («)-$ (M) + (l~Pi)$ (u) 

J ° V y/l=tf 

— <[>( ii )+(1 — pi )$“(«) 


dv 

dv 


+ /« “ $ ( 

u ~k- I-oc $ (^+^=f $_ (w))</'(^V 1 'P? + $_ ( u )) yj^-pfdv 

+r(i-K 


v / 1 - Pi + $ (u)W 1 - p? dv 


V 1 - Pi (/-oo $ W( $ ~ (w)) dv + 0(V1 - Pi)) 
+ V /l -Pi(/o °°( 1 - $(E)<K $_ (ii)) dv + 0(V1 - Pi)) 

= ^ i ( vk<A( $ “( u )) + °(Vi - p*)) 

+V 1 - Pi (vb^ - ^)) + °(V 1 - ft)) 

= u— \ + 0(1/log ra) 

fo(v - ±) 2 c(u,v; Pi) dv 
= \ - 2/ 0 1 (v - i)Oi(u,v;pi)dv 


JO V 2J 

= (v - |) 2 - 2 / 0 > - |)$( ri,H y |rM ) dv 


+ 2 />-^ 


l)(i-*( 


2 ' v 

<J> _ (v)— <I> - (li) + (l— pi)<E> - (w) 


v/w 


dv 


= (« - \Y - 2 V 1 -P 1 (/_ oo( v ~ §)$(«)<£($“ («)) dv + 0 (y/l - pi)) 

+ 2 \/l — pf ( /,T( W - 5 H 1 ~ $(v)) 0 ($~(u)) dv + 0(^1 - Pi)) 

= 0 - |) 2 - 2 \A - pf + i)^( $_ ( u )) + °(v /1 - ft)) 

+V 1 -P?((3 - ^)/>($”(m)) + 0(^1 - pO) 

= (w - |) 2 + - P 2 i + 0 ( 1 /log n). 


By (5.10), (5.11), C(u, v; 1) = u A v and 

—C(u,v;p) = - , 

rfp 2^v/l^ 


; exp - 


(4> (it )) 2 — 2 p$ (m)< 1 > (v) + (<h (v)) s 


2(1 -P 2 ) 


(see Plackett (1954)), we have 


E%=E% 

= Jrdtw E"=l ELl In fo fo fo ( U 1 A U 2 - «1«2)(«1 - |)(«2 ~ f)* 
c(u lt v 1 ;pj)c(u 2 , v 2 ; pu) dvidv 2 du 1 du 2 
= fg /„(«! A u 2 - viu 2 )(mi - |)(it 2 - I) duidu 2 + 0(1/log ra) 

= ^ + 0 ( 1 /log n), 


(5.10) 


(5.11) 


(5.12) 


(5.13) 
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E% 


= fo fo( u ~ 5 ) 2 c{u,v,Pi) dvdu - (iarcsin(f)) 

= fo( u ~ D 2 (( u - s ) 2 + </>($”(u))\/l - Pi) du + 0 ( 1 /log n) - (^arcsin(i)) 

= §S + r g n /K) fo( u ~ \) 2 <t>{<S>~{u))du- yij +0(1/logn) 

_ j_ V 2 y/m(i/n) fl( u _ l)2^($-( u ))du + 0( 1 /logn), 


— -L 

180 ^ 


E(Zi,iZi, 2 ) 

= (^IF E?=l ELi fo fo fo fo (C( u l, v 2; Pi) - U lV2 )( Vl i)(U2 - |)x 
c(zti, «i; pj)c(u 2 , v 2l /+) dv\du 2 du\dv 2 
= fo 1 fo( Ul Av 2 - Uiv 2 ){ui - \){v 2 - i) duidv 2 + 0 ( 1 /log n) 


720 


+ 0(1/log 71) 


and 


E(2'i i i^j ) 3) — E(Z ii2 Z ii3 ) 

= 0 +T E"=l fo 1 fo 1 fo 1 fo 1 ( f ( U 2 ^ U l) - M l)(^l - |)( U 2 - s)(«2 - §)x 

c(iti, Mi; pj)c(u 2 , m 2 ; Pi) dv\dv 2 du\du 2 


= f 0 fo (I{u 2 < Ml) - mi)(m 2 - |)(mi - |)(m 2 - 7 j)dMidM 2 + 0(1/logn) 
= “300 + 0 (i/logn). 

Hence, it follows from ( 5.13[ )-( [5.16[ ) that 

1 Fez 2 = ^fo^+~^ ds f 1 

71 1 


i=l 




(/>($ (m)) dM + 0(1/ log ti). 


It is easy to check that 

E (l E(^ 2 -E^ 2 )] = ^E(l 2 ^El 2 ) 2 = 0(l/n), 

\ 2=1 / 2=1 

which, combining with ( 5.17| ), implies that 

4 W+^*) (/'<“- I)V(*-W)*.) 


2=1 

Obviously we have 


max 

l< 2 <n 


(log Tl ) 1 / 4 




A 0 and E ( max 


(log 77) 1 / 2 


yl< 2 <n Tl 


Z? =o(l). 


(5.14) 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


(5.19) 


Hence, it follows from ( |5.9[ ), (5.181, (5.19) and Theorem 3.2 of Hall and Heyde (1980) that 

—Zni(<*,/3,7) -» N ^ 0 , yja + /3s^ ds^j ( J (u - ^) 2 (j>($~(u)) dvfj^J . (5.20) 

Note that the above limit has a nonstandard rate, which can be explained as follows. When pi = 1 for 
i = 1, ■ • • , 77 , we have 


u«,/ 3,7) = E ((to) - ^) 2 - ^) = E (( 


- - -) 2 - — ) 

71 + 1 2 2 12/’ 
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which becomes a constant. However, l n2 (a, 0,j) and ^3(0:,/3,7) are non-degenerate due to the involved 
factors ( i/n ) 7 and ( i/n ) 7 log(i/n). That is, deriving the asymptotic limit of l„ i(a, /3,7) needs finer expansions 
than the other two quantities. Below we show the asymptotic limits for both l n2 and l n 3 have the standard 
rate l/yfn. 


Define 

z* = p^T) Ej=i(n) 7 fo ^ u ) ~ u )( v - \) c ( u , v; pj) dudv 

+ (^T) Ej=i(n) 7 fo < v)-v){u- \)c{u,v,pj)dudv 

+ {Z 1 - ^ arcsin (f))(^) 7 
= Z* 1 + Z* 2 + Z* 3 

1 

= pi+i) Ej=i(^) 71o g(n)/o 1 /o 1 ( / ( [/ i <u)-u)(u- ^)c(u,v; Pj) dudv 
+ (+FT) E"=i( a ) 7 log(^) Jo (I(Vi <v)-v)(u- i)c(w, v,pj) dudv 
+ { z i - ^arcsin(f ))(^log0 

_ ^** _j_ —j— * 


Similar to the proof of (5.91, we can show that 


f ^4„ 2 (a, /3, 7 ) = +°p (VV") 

1 z n3 (a,ft 7 ) = -4 ELi^r + o P (i/Vn). 


Like the proofs of (5.13) (5.16), we can show that 


(5.21) 


E ' 7*2 _ XT' ' 7*2 

^i,l — ^ 2 

= (^Ei=i(^) 7 ) 2 /o L /o 1 ( u i Au 2 -wiw 2 )(wi- |)(u 2 - i) du\du 2 + 0(1/logn) 

= 720(1+7)^ +Q( 1 / 1 °g^) 1 


I EZ*i = ^(1)27 + 0(1/^10^), E(Z* 1 %) = ^ 0( ^ F +O(l/logn) 

I E(%Z* 3 ) = E(Z* 2 Z* 3 ) = + O(VVIogn), 


jg ^**2 _ g ^**2 

= (iT Ej=i(n) 71 °g(^)) Jo fo(uiAu 2 -u 1 u 2 )(ui - I )(u 2 - §) du\du 2 + 0(1/ logn) 
= 720(1+7)* +0(1/logn) 

and 

f E%3 2 = T | 5 ( A )^log 2 (A) + 0(l/Vb^), E(%Z?S) =™(T+^ + 0(1/logn) 
1 E(%Z** 3 ) = E(Z*; 2 Z%) = 3 W ^ f (A) 7 i og (i) + 0(1/Vb^), 

which imply that 


E( 


1 sr^ n r?*\2 _ 1 7* 

y/E 2^i=l J ~ n 2^i=l 1 

7**\2 _ 1 


?* 2 + zi 3 ) 


i,3J 

■** \2 


e (^ e”=i ^**) 2 = J Eh m:;i + z** 2 + z**) 


180 ( 1 + 27 ) 

, _1_ 

90(l+27) 3 


180(l+7) 2 

_1_ 

" 180(l+7) 4 ’ 


(5.22) 


Like the proof of (5.20), by using (5.22), we can show that 

-^l n2 (a,/3,j) -+ N fo, 


180(1 + 2 7 ) 180(l + 7) 2 
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and 


-^/ n3 (a,/3,7) 4 7V(0 

Some further tedious calculations show that 

(log n ) 1/4 


’90(l + 2 7 ) 3 180(1 + 7 ) 4 


and 


EIL i E (W) 


(logn ) 1 / 4 


(log n ) 1/4 


+ - 


E n V ** 3 Tp/y \ 

i=1 Z^j=1 Z^/c=l 


E?=i ( 


+ 


=1 1 720(1+7) ^ 720(1+7) 360 V n 


_3-('■1)7 

360 V n > 


'720(1+7) 

O(l), 


+ 


720(1+7) 


\l /4 


_— ( — )1 - 

360 'n' 


360(1+7) 360(1+7) “ 180 U 


+ t+(+ 7 + 0(l/\/logn 


ik ^EIUE(W) 

E n \ ^ 3 \ +3 ( 'y ' 7 ** \ 

2—1 2-^=1 2 ^ fc=l 


(log? 


xl/4 


Ei=l ( 720(l+7) 2 720(l+7) 2 360 ( n ) 7 l°g( „ ) 

— 720(l+7) 2 _ 720(l+7) 2 ~ 360 («) 7 1°§( + + 360(l+7) 2 + 360(l+7) 2 

+t§o+) 7 lo s+) + o(VVEg^)) 

= 0 ( 1 ) 

sELiE (Z*z**) 


_ 1 \“^ n \ ^3 \ ^3 t7'/ > 7 * 7 * 

— n Z^2=l 2-^k=l 


* l ) 


B Et=l ( 720(1+7)3 720(1+7)3 360(1+7) (n ) 7 l°g( n ) 

_I_1_I_('ll l OB +i\ -1_1_(1)7 _I_("1)1 

720(1+7)3 720(1+7)3 360(1+7) t n > n > ~ 360(l+7) 2 1 n > T 360(l+7) 2 'nl 

+ TM(^) 271 og(++0(l/Vbg^)) 

180 (l+ 27) 2 + °( 1 )- 


180(1+7)3 

Hence, by Cramer device, we can show that 


/ (log n) 4 / 4 1 1 

- T= —7), —F=ln2{oi,p,i), — r l n3 (a,P,'y) 

V V n v n v n 


iV(0,E). 


(5.23) 


It is straightforward to check that 


log n dl n i(ot)f3tf) . \/3 log n dl n i(a,^,j) , \/3 

n 0 a 67 t ’ n 0/3 67 r(l+ 7 ) ’ 

log n 0^ 7 i 1 (a,/3, 7 ) __ -y/3/3 log n 0 Z ti 2 (Q!,/3, 7 ) _v \/3 

n O 7 67 t( 1 + 7) 2 5 n 0a 677 ( 1 + 7 ) ’ 

logn dl n 2 (a,P, 7 ) v logn 0 ^ 2(03 P, 7 ) v_ \/3/3 

n 0/3 677 ( 1 + 27 ) 5 n O 7 677 (l+ 27) 2 ’ 

logn dl n 3 (q,/ 3 , 7 ) -y/3 logn 0 ^ 3 ( 0 ,/3, 7 ) _ \/3 

n 0a 677 (l+ 7) 2 5 n 0/3 677 (l+ 27) 2 ’ 

logn 01 ^ 3 (a,/3, 7 ) y/3/3 

n O 7 377 (l+ 27) 3 ' 


Hence, the theorem follows from (5.23), (5.24) and Taylor expansion. 


(5.24) 


□ 


Proof of Theorem [3J As in the proof of Theorem [2] we define 


Zi 


= J^»T.USofim%WU,<n) 
+ph)E?„/„7o^, (HVi<v) 
+ (•!'-(V,) - p, 

= Zi t 1 + 7,72 + + 3 j 


u)c(u, v; Pj) dudv 
- v)c(u , u; pj) dudv 
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z* = (^xy E "=1 (V 7 fo fo *(*-(«)) ( J ( Ui ^ u ) - u)c{u,v, Pj )dudv 
+^TI) E "=1 (i) 7 fo fo ***% ( J ( y » ^ «)-v)c(u,t;;pj)d«^ 


+ 


Z * _l y* _i_ y* 

i,l ' ^,2 ' ^i ,3 


and 


^ 7 ** _ 


Since 


we have 


= (STI) E"=i(^) 71 °g(^)/ 0 1 fo' +**-$)) ( f ( u i ^ u)-u)c(u,v;pj)dudv 
+ (STT) Ej=i(£) 7 log(i) /q 1 fo 4 **-%j (I(Vi <v)~ v)c(u,v, Pj ) dudv 

+ ^-(U l )^-(V z )-p l )(fr log(i) 

_ _l_ _(_ 

/ ^ _ (i;)c(w, v; dv = p* $“(?/) and / <£>“ ('it)c('W, v; p*) du = 

Vo Vo 


^,1 = ( dTT E"=i Pi) fo' «*-$)) ( f ( U i ^ u ) - w ) dM - 
^,2 = f E"=l Pi) fo 1 #(*-$)) ( f ( V i < V )~ V ) dv > 

K i = (dr e;=i Pidr) fo 1 <«) -«) d«, 

% = (dTT E” =1 Pi(^) 7 ) fo 1 £tyj(W <v)-v) dv > 

Z*f = ( E "=1 Pj(n) 7 lo g(„)) fo 1 #*-$)) Wi < u )^ u ) du, 

Z% = (dn E"=i Pi(£) 7 log(i)) fo 1 tP$jj(W <v)~ v) dv. 

It is straightforward to check that 

E % = (nTT E"=l Pi) fo 1 fo 1 (£))?($-(u a )) (m A U 2 - « 1 « 2 ) 

= (srr E"=1 Pi) 2 2 fo 1 f 0 Ul ^-(l" 1 )j*(J-^ 2 )) M2(l - Ml) du 2 dui 
= (dTT E"=l Pi) fo 1 M*-(ul)) (! - M l) f-J Ul) 

= (nil E"=l Pi) fo 1 »(*-fri)) ^ ~ Ul ){( $ ~( Ul )) 2ui + /-oo (Ml) M 2 #M} dm 
= (^E"=l Pi) fo 1 MZ-iulv i 1 -Ml){(^~(»l)) 2 Ml + $~K)<K$'~(Ml)) -Ml}dM! 
= (dTT E"=i Pi) fZo “(! - $(m)){m 2 $(m) + u0(u) - $(u)} du 
= - (sti E"=1 Pi) JZo u ( X “ $ ( u )) 

= (dTT Ei=i Pi) fZo < ^( u )i 1 - $ ( M ) - W>(m)} du 
= I (s+lE"=i Pi " 

by noting that u(l — $(u))$(u), u 3 (l — $(u))$(u) and u(j> 2 {u) are odd functions 

\ 2 


1 


1 


e ^ 2 2 = d I I and EZ? 3 = l + 2p*-p* = l + p*. 


n + 1 


i=i 


By (5.12), we have 


$ (u)$ (u) dC(u,v;p) 
o A) </>($” («))</>($” (m)) dp 


dudv = p 


(5.25) 


(5.26) 
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for any p £ (—1,1). Taking derivative with respect to p at both sides, we have 


$ (u)$ (v) d 2 C(u,v; p) 
o Jo </>($" (u))<l>(®-(v)) dp 2 


dudv = 1 


for any p £ (—1,1). Therefore 


= So fo ( u A v-uv) dudv 

= fo fo („)) ( C ( u > v ’ X ) - uv ) dudv 

= fo 1 fo 1 ft*-(%)*(*%)) ( C ( u > v > - uv ) dudv 

+pi( 1 - Pi) + |(1 - Pi) 2 +o(l/log 2 n), 


which gives 


E( + , 1+2) 


Since 


( 1 

S+T 2-<j=i Pi 
E(I(Ui <w)$-(C/i)$-(Fi)) 


= (hrr E"=i Pi) ^ fo 1 fo 1 w-ffwJ-lv)) (°( u ’ v ; ft) - uv ) 

= (s;iE"=ift) (5 -Pi(l-Pi) - 5(1 -P*) 2 ) + o(l/log 2 n) 
+ o(l/log 2 n). 

= E (/(C/i < u)$-(17 i )E($-(yi)|$-(C/i))) 
= E(I(U i <u)$-(U i )p i $-(Ui)) 

= Pi f*J U) v 2 0(v) dv 
= -Pif-J U) vd(l>(v ) 

= “Pi ($~ («)</>($"(«)) - u) , 


we have 


E(Zi,iZi,3) — E(+2+3) 


(fAi E "=1 Pi) fo 1 0(V ( pV * (“$ («)) + M - w) du 


- (F_ 

\n+l 


E”=! Pi) Pi- 


Put p = n 1 Ej=i Pi • Then 

I P = 1 ~ a+ l / g ( n +7) + oCV^g 2 + 

l k E?=i P 2 = ^ E"=i(Pi - P) 2 + P 2 = i^(i^ ^ (T^f) + P 2 + o(l/log 2 n). 
It follows from ( 5.25| H5.29) that 

A E"=i E + 2 

= (sir E "=1 p,) 2 + 1 + ^ E ?=1 pf + £ EIU p 2 (^t E "=i p,) 2 

-£ E?=1 Pi(iiTT E"=1 Pi) + 0(1/log 2 n) 

= p 2 + 1 + p 2 + i E"=i (Pi - p-) 2 + (P“ 2 + £ E ”=1 (Pi - p-) 2 )p' 2 

—4p 2 + o(l/log 2 n) 

= 4 ( Q+ fog ( n +7) ) 2 + iIp+IaW “ ( 1 + 7 ) 2 ) + o(l/\og 2 n). 

Similarly we can show that 

^E”=iE^ 2 

_ 1 V^ n \"^3 Y^3 Ti'fV* C7* 

— „ 2^i=l+j = l 2^fc=l 

— - V" ( 1 „ .li¬ 

ra 2^i=l y 2(l+7) 2 ^ 

1 1 _ri__, 

1 + 7 ' 


r 2 (l + 7 ) 2 ^ 2 (l+ 7) 2 I +7 

2 2 ' 0(1), 


1 v . _l_f'A'|7 

2 (l+ 7) 2 I +7 \n> 

_1_('A'\7__l_('i’i7 4 . of A 

1+7 in/ I+ 7 in/ 1 + 7 'n/ ' 'n 


2 ( 1 ) 


I +27 (I+ 7) 2 


(5.27) 


(5.28) 


(5.29) 


(5.30) 


(5.31) 
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1 


— E Z**^ 

_ 1 V^ n T 7 ( 'y** ' 7 ** \ 

— n Li=l Lj=l Lfc=l 

= sElU ( 2(1771 + 2 (Tq 7 F + (i+TF(n) 71 °s(s) 

+ 2(1+77 + 2(1777 + ( 1 777 (n) 7 lo S(ra) 

+ ( 1 + 1 7 ) 2 (n) 71 °g(n) + (1777 (n) 7 l0g (n) + 2 ( S lo § 2 ( ^ )) + °( 1 ) 

_ 4 2 


(l+2 7 )3 (I+ 7 ) 

logn 

2=1 


°( 1 ), 


E E (^n = 0(1), 1 °^ n E E &%*) = °w 


n 


and 


sE^W) 

1 V'' n V^3 17 / ry* ry** \ 

n Z^i =1 1 Z^fc=1 

n E)i= 1 ( 2 ( 1+77 “ 2(1777 “ lT7^n ) 7 ^ 0 g ^n) 


" 2 ( 1++) 3 2 (l+ 7) 3 1 +-V (n ) 7 ^°S^r 7 , ) 


+ 


.(A)7 + 
V n / 1 


(l+7) : 

2 2 

(1+277 + 0+7? 


1+7 ' n 

1 (7) 7 +2++log+ 


(i+7) 


■o(l) 


o 


(!)• 


Therefore, using (5.301 (5.34) and the same arguments in proving (5.23), we can show that 

( ^0(«» ft 7 ), 7 ^ 2(7 ft 7), -^ 3 +, ft 7 )) 4 iV(0, E*). 

V v n v 71 v n / 

It is straightforward to check that 

logn 9(*i(a,/3,7) 1 logn <91* 1 ( CK >/+7) 1 logn 51*!(a,7, 7 ) 7 

n 5a 1 n 57 1 + 7 ’ n 57 (I+ 7) 2 5 

logn 51* 2 (a,/3, 7 ) 1 logn dl* 2 (a,/3, 7 ) . 1 logn 51* 2 (a,7, 7 ) , _ 7 

n 5a 1+7 ’ n 57 I+ 27 7 n 5+ (I+ 27) 2 ’ 

logn 5 t* 3 (a,/ 3 , 7 ) 1 logn 51* 3 (a,ft, 7 ) 1 logn dl* 3 (a,/3, 7 ) _ 27 

n 5a 


(I+ 7 ) 2 ’ n 


57 


(I+ 27) 2 > n 


57 


(I+ 27) 3 ' 


Hence, the theorem follows from (5.35), (5.36) and Taylor expansions. 


Proof of Theorem [7} It follows from the proof of Theorem [2] with known 7 = 1. 

Proof of Theorem [5| It follows from the proof of Theorem [3] with known 7 = 1. 

Proof of Theorem [6| Note that 

(log n)Q”(s) -+ — v ^ m ( s ) anc j cos(27 tQ(s)) -+ — 
67t 2 


and 


E ((7i(Xi) - |)(F 2 (Fi) - |) - + arcsin(f)) 

= E ((i+X,) - i) 2 + 2 (17) - l) 2 ) - (£ arcsin(f )) S 

-+ 1 — (JL) 2 = 

^ 80 1121 


1 

180' 


(5.32) 


(5.33) 


(5.34) 


(5.35) 


(5.36) 

□ 

□ 

□ 

(5.37) 

(5.38) 


It follows from (5.38) and the standard arguments in local linear estimation (e.g., Fan and Gijbels 0) that 
Vnh ^Q(a) - Q(s) - 1 Q"{s)h 2 J t 2 k(t) dtj -%n(o,^J k 2 {i) dtj . (5.39) 


20 













































Hence it follows from (5.37) and (5.39) that 


^ (m(s)-m(s)) 

= — 4n cos(2TrQ(s))Vnh(Q(s) — Q(s)) + o p ( 1) 

-4- N m"(s) t 2 k{t) dt, k 2 (t ) dtj , 

i.e., the theorem holds. □ 

Proof of Theorem [?| It follows from standard arguments in local linear estimation. □ 
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